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Meditation proficiency is related to trait-like (learned) effects on brain function, developed 
over time. Previous studies show increases in EEG power in lower frequency bands 
(theta, alpha) in experienced meditators in both meditation states and baseline conditions. 
Higher gamma band power has been found in advanced Buddhist meditators, yet it is not 
known if this occurs in Yoga meditation practices. This study used eLORETA to compare 
differences in cortical source activity underlying scalp EEG from intermediate (mean 
experience 4 years) and advanced (mean experience 30 years) Australian meditators from 
the Satyananda Yoga tradition during a body-steadiness meditation, mantra meditation, 
and non-meditation mental calculation condition. Intermediate Yoga meditators showed 
greater source activity in low frequencies (particularly theta and alphaD during mental 
calculation, body-steadiness and mantra meditation. A similar spatial pattern of significant 
differences was found in all conditions but the number of significant voxels was double 
during body-steadiness and mantra meditation than in the non-meditation (calculation) 
condition. These differences were greatest in right (R) superior frontal and R precentral gyri 
and extended back to include the R parietal and occipital lobes. Advanced Yoga meditators 
showed greater activity in high frequencies (beta and especially gamma) in all conditions 
but greatly expanded during meditation practice. Across all conditions (meditation and 
non-meditation) differences were greatest in the same regions: R insula, R inferior 
frontal gyrus and R anterior temporal lobe. Distinct R core networks were identified in 
alpha 1 (8-10 Hz) and gamma (25-42 Hz) bands, respectively. The voxels recruited to these 
networks greatly expanded during meditation practice to include homologous regions of 
the left hemisphere. Functional interpretation parallels traditionally described stages of 
development in Yoga proficiency. 
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INTRODUCTION 

The introduction of Yoga and Buddhist meditation practices into 
Western countries has extended their application beyond their 
original spiritual goals (Shapiro, 2003) to include a range of 
mind-body interventions for health-related problems (Ospina 
et al., 2007). Research has explored the neurophysiology of 
meditation during the actual practice (as examples of specific 
altered states of consciousness) and changes persisting into non- 
meditation conditions (as examples of neuroplasticity). A central 
question in this investigation is how the level of meditator profi- 
ciency contributes to the development of these effects. 

As yet there is no accepted objective measure of meditator 
proficiency. Most studies report proficiency in years, although a 
more recent trend uses the more accurately calculated "hours of 
practice," "based on daily practice and time spent in meditative 
retreats" (Brefczynski-Lewis et al., 2007). The term "advanced" 
has usually been reserved for meditators with more than 20 
years' experience (Arambula et al, 2001). Numerous studies 



conducted over the last 40 years with Western meditators, usu- 
ally with less than 10 years' experience, have reported increased 
power and coherence in the alpha and theta frequency bands 
during meditation practice (Cahn and Polich, 2006). Striking 
increases in gamma band power have also been reported in stud- 
ies with "advanced" Buddhist monks as meditators (Lehmann 
et al., 2001; Lutz et al., 2004). A recent study also found sig- 
nificantly increased gamma power at parieto-occipital electrodes 
in a group of "advanced" Western Vipassana meditators (with 
mean 20 years' meditation experience — although the nature of 
that experience is undefined) engaged in a "mindfulness" body- 
scan meditation when compared to a deliberate mind-wandering 
instruction condition (Cahn et al., 2010). 

Fell and colleagues proposed that meditators from different 
traditions progress through similar developmental stages, which 
are marked by changing EEG frequency patterns (Fell et al., 
2010). Initial expertise is reflected in changes in slower (specif- 
ically theta and alpha) frequency bands, but these effects are 
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not considered to be meditation-specific. The authors propose 
that an "advanced" stage, only reached by experts, is marked 
by increased synchronized (fast) gamma band activity, related to 
"processes of cortical restructuring and learning" which facilitate 
"specific meditation-related states of consciousness," with unique 
electrophysiological signatures. 

A recent focus in EEG and brain imaging meditation studies 
has been on the pivotal role of the modulation of the "default 
mode network" (DMN) in the understanding of meditation 
training effects (Brewer et al., 201 1; Berkovich-Ohana et al., 2012; 
Malinowski, 2013). This neural network is considered to be active 
when effortful attention is not required to direct responses to 
the external environment. Anatomically the DMN consists of two 
major interacting subsystems defined by the cluster of activa- 
tions and functional connections around medial prefrontal and 
medial parietal hubs, respectively (Buckner et al., 2008). A series 
of both EEG and fMRI studies have converged on the finding of 
reduced activity in these major hubs of the DMN during Buddhist 
and/or mindfulness meditation practice (Farb et al., 2007; Holzel 
et al, 2007; Berkovich-Ohana et al, 2012). The DMN is known 
to be activated by spontaneous self-related mentation, directed 
to either the future or the past (Schooler et al., 2011), character- 
istic of distraction or mind wandering. Both Yoga and Buddhist 
meditation traditions 1 note the emergence and progressive reso- 
lution of distraction and mind wandering during the early phases 
of mediation training irrespective of content. Rather than belong- 
ing to the essence of states cultivated by respective meditation 
practices, current DMN findings may only reveal the expectable 
psychological effects of early steps in the process of learning to 
meditate. 

Recent brain imaging techniques have begun to reveal train- 
ing effects of a variety of meditations, across a range of practice 
periods, on the neuroplastic response of the brain in both gray 
matter and white matter structures (Lazar et al., 2005; Pagnoni 
and Cekic, 2007; Holzel et al, 2008, 2011; Luders et al, 2009, 
201 1, 2012a,b; Vestergaard-Poulsen et al., 2009; Grant et al., 2010; 
Tang et al, 2010; Murakami et al, 2012; Kang et al., 2013; Leung 
et al., 2013). In EEG studies Berkovich-Ohana et al. (2012) found 
significantly greater gamma-band activity at rest in mindfulness 
meditators than controls in a cluster of electrodes in the right 
parieto-occipital region and the reverse in a cluster of electrodes 
in the right frontal region, however no significant relationship was 
reported for these measures with meditation experience within 
the mindfulness meditation group. Cahn et al. (2010) found 
gamma power at occipital electrodes increased significantly from 
rest to "vipassana meditation" in those with greater than 10 years 
of daily practice but only marginally in those with less than 10 
years of daily practice. Although their findings vary widely (as 
might be expected from their methodological differences) these 
studies point to the importance of evaluating both the long-term 
impact of specific meditation practices and the way in which 
that impact unfolds over time. Understanding these medium and 



In the present study Buddhist and Yoga meditation practices are neither 
assumed to be equivalent nor different at a neurophysiological and/or psy- 
chological level, rather this is something to be established by empirical 
investigation. 



longer term changes may require a focus beyond the role of 
the DMN. 

Within the Yoga tradition sequential stages in the develop- 
ment of meditation proficiency are delineated in the "eight-fold" 
path of Patanjali's Yoga Sutras (Radhakrishnan, 1999). Following 
four "external" stages, the practitioner progresses through the 
"internal" stages of "pratyahara" (sense withdrawal), "dharana" 
(concentration), "dhyana" (absorption), and finally "samadhi" 
(self-realization). Repeated experience of these states during med- 
itation practice leads to a long-term progressive refinement in 
the "sense of self" (asmita). One's self- identity becomes pro- 
gressively detached from identification with externally-oriented 
perceptions, then from identification with the body, then from 
identification with thoughts, to absorption in the object of med- 
itation in samadhi. The process has been likened to peeling 
the skins from an onion (Maheshwarananda, 2005). However, 
the electrophysiological changes corresponding to such a devel- 
opmental sequence will not be evident within current cross- 
sectional designs comparing only short-term and non-meditators 
or short-term and long-term meditators. While longitudinal 
studies are ideal, as a first step it is necessary to track the transition 
between medium term and long term effects of Yoga meditation 
training. 

The present study made use of the long period of estab- 
lishment of Satyananda Yoga in Australia for the availability 
of "advanced" Western practitioners with over 30 years' expe- 
rience. The study compared "advanced" Australian Satyananda 
Yoga teachers (SYT) with students studying to become SYT 
(having an intermediate level of experience). Two meditation 
practices were used for the study. The first "kaya sthairyani" 
is a preparatory body-steadiness practice designed to take the 
practitioner into the "pratyahara" stage (withdrawal of the mind 
from the external world). The second practice "japa" uses men- 
tal repetition of a personal mantra to move from "pratya- 
hara" to the "dharana" stage (a focussed internal awareness 
of the mantra). Advanced practitioners may progress further 
to "dhyana" (absorption in the mantra) or even "samadhi" 
in this practice (Saraswati, 1983). Based on this tradition, we 
hypothesized that student meditators are more likely to expe- 
rience sense withdrawal; the advanced practitioners are more 
likely to progress to the later stages when engaged in these 
practices. 

In keeping with results from moderately experienced Western 
meditators and the demands of sensory inhibition, we hypothe- 
sized firstly that the Satyananda students (with an intermediate 
level of meditation practice) would show more EEG activity in 
the lower frequency bands (theta, alphal, and alpha2). It was 
expected that this activity would be higher in the students than 
in the SYT group as a consequence of the teachers' development 
beyond this stage of practice. Secondly, it was expected that the 
SYT would show greater activity in the higher frequency bands 
(beta and gamma) compared to the Satyananda Yoga student 
group in both meditation conditions. Thirdly, extended prac- 
tice (spanning years) in both student and teacher groups was 
expected to result in enhanced neural connectivity and thus trait 
activation in networks habitually activated in each group during 
meditation outside the context of meditation itself. Fourth and 
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finally, to the extent that Yoga meditation practices engender sim- 
ilar DMN changes to the Buddhist-mindfulness practices recently 
investigated by meditation researchers, it was expected that the 
cortical foci of these length of practice differences would lie in the 
anterior and posterior midline hubs of the DMN at least during 
active meditation practice 2 . 

MATERIALS AND METHODS 
PARTICIPANTS 

Twelve Satyananda Yoga practitioners (teachers and students) 
were recruited for the study from the Satyananda Yoga Academy 
in Mangrove Creek, NSW, Australia by the first author. At the time 
of the study, he was a resident lecturer at the Academy. In his 
sixties, he is a karma sannyasin disciple of Swam Satyananda of 
Munger, India. Participants were recruited initially by direct con- 
tact for a qualitative study of meditation. Following that study, 
participants for the EEG study were selected from suitable volun- 
teers. All the teachers had received initiation as sannyasin disciples 
of Swami Satyananda and had been regular meditators for over 
20 years. The student group was studying an accredited course to 
qualify as SYT All participants had received a personal mantra 
from their guru (Swami Satyananda or his successor, Swami 
Niranjanananda), which they used in this study. In pre-study 
interviews, all reported they were free from medical, psychiatric, 
or drug usage issues that might alter their brain functioning. 

The participants were divided into two groups — SYT, three 
male and three female, age range from 44 to 63 years (mean = 
54 years, SD = 6.5 years) and Satyananda Yoga students (SYS), 
three male and three female trainee Yoga teachers, age range 30- 
51 years (mean = 42 years, SD = 8.0 years). The teacher group 
was significantly older [f(nj = 2.90, one-tailed, p < 0.05]. The 
SYT group had a mean of 30 years regular practice (range: 24-37 
years) and the SYS group a mean of 4 years (range: 3-5 years). 
Based on an average of 1 h regular practice a day, this would 
equate to a mean of 1 1,000 h for the SYT group and 1500 h for 
the SYS group. 



2 The research reported here formed part of the Ph.D. thesis of the first author, 
supervised by the second and third authors. Following the framework of Cahn 
and Polich (2006) the thesis sought to determine expected trait (length of 
practice) effects, state (meditation vs. non-meditation) effects and state by 
trait interaction effects of Yoga meditation practice in EEG cortical source 
activity. Specific trait expectations (see above) were principally based upon 
the work of Farb et al. (2007) (for cortical regions) and Fell et al. (2010) (for 
frequency bands) and it is the results pertaining to these hypotheses that are 
reported here. It was expected that state effects would be found in higher high 
(i.e. gamma) frequency source activity in cortical regions specialized in pro- 
cessing the different contents of specific meditations (e.g., bilateral parietal 
and somatosensory cortex for body scan and left inferior frontal gyrus for 
mantra) during those specific meditations as compared to the non-meditation 
condition. Following the person by situation interactionist approach in per- 
sonality psychology it was expected that state by trait interaction (as distinct 
from additive) effects would be highly significant (and informative) although 
the specific form of those effects was not predicted. Neither the expected 
state nor state by trait interaction effects were found. We consider these null 
findings to be highly informative (although beyond the scope of the current 
report) and they have reorientated our perspective to focus on the under- 
standing of trait changes arising over the long term time course of meditation 
practice. 



SETTING 

The study was conducted in a small meditation room at the 
Satyananda Yoga Academy to provide an "ecologically valid" sit- 
uation conducive to the attainment of deep meditation states. 
After being fitted with the Compumedics 32 channel EEG "Quik 
Cap," participants sat in their usual cross-legged meditation posi- 
tion on the floor, supported by cushions, with the room dimly 
lit by a candle on a meditation table. The small battery-operated 
Compumedics "Siesta 802" recording unit attached to the "Quik 
Cap" was the only electronic device in the meditation room. This 
unit transmitted signals by radio to a laptop computer in an 
adjacent room. 

PROCEDURE 

To obtain the most authentic meditation experience possible, the 
sequence of conditions was selected to resemble the practitioners' 
usual meditation practices. We considered that a design incor- 
porating a counterbalanced order of conditions would introduce 
a conflict with traditional practice, as kaya stairyam would usu- 
ally precede but not follow japa (Saraswati, 1981). Following the 
mental calculation condition, the participants performed four 
meditation practices in the same order (for brevity only the first 
two of these are reported here). 

The experimental conditions were: 

• Non-meditation condition — "Calculation" — mentally count- 
ing backwards from 200 by 4 s — (5 min). 

• Meditation 1 — "Body-steadiness" — Satyananda Yoga practice 
"kaya stairyam" — awareness focused on body-steadiness and 
awareness of flow of the natural breath — (5 min). 

• Meditation 2 — "Mantra" — Satyananda Yoga practice "japa" — 
mental repetition of personal mantra, using mala (beads). 
The mantra consisted of a short Sanskrit phrase — (10 min) 
(Saraswati, 1981). 

All conditions were performed, sitting in cross-legged medita- 
tion posture, with eyes closed. Self- report ratings of "meditation 
depth" for each condition were obtained immediately following 
the EEG recording. These ratings were made on a ten-point visual 
analog scale, based on the "Meditation Depth Questionnaire" of 
Ulrich Ott (2001). The zero point was "no meditation state" and 
ten was "deepest meditation I have experienced." 

The study was approved by the RMIT University Human 
Research Ethics Committee. 

EEG DATA COLLECTION 

EEG signals were obtained using a Compumedics "Quik Cap" 
from 25 scalp electrodes, based on the International 10/20 system, 
referenced to left mastoid. Electrodes were placed at FP1, FP2, F3, 
F4, C3, C4, P3, P4, Ol, 02, F7, F8, T3, T4, T5, T6, Fz, Cz, Pz, Oz, 
TP7, TP8, CP3, CP4, FC3. Four additional channels were allo- 
cated to eye movement detection, with electrodes positioned on 
the outer canthi of each eye and above and below the left eye. 
The sampling rate was 256 Hz. Data was acquired via radio sig- 
nal to a laptop computer running Compumedics Profusion EEG 
software. 
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EEG DATA PRE-PROCESSING 

Following the recording session, the data was exported from 
the Profusion EEG in EDF format for input into the EEGlab 
(Delorme and Makeig, 2004) program running in Matlab. Data 
was preprocessed through the FASTER algorithm (Nolan et al., 
2010) which employs independent components analysis (ICA) to 
identify and remove both physiological (eye movements/blinks, 
muscle movement, skin potentials) and non-physiological (elec- 
tromagnetic interference, electrode pop offs and drift, shifting 
electrodes and residual white noise) sources of artifact from the 
recorded EEG. FASTER interpolated missing or bad channels, re- 
referenced to the common average and applied a bandpass filter 
of 1-45 Hz to remove drift and further high frequency artifacts. 
FASTER detects and removes ICA components with properties 
uncharacteristic of cortical signals or conversely with properties 
characteristic of specific artifact sources. Z-score thresholds for 
rejecting artifactual components were set at 3.0 (except for eye 
movement, which was set at a threshold of 1.8). EEGLAB version 
9 (Delorme and Makeig, 2004) was employed for criterion-based 
artifact rejection of epochs with values greater than ±75 mV. The 
recorded EEG was then subject to visual inspection as a final check 
of artifact removal. 

Ten sample epochs, each of 5 s duration, were extracted for the 
first 5 min of each condition for each participant from the cleaned 
datasets, commencing 100 s from the start of the practice and 
then at 20 s intervals. If a selected epoch showed residual artifacts 
on visual inspection, the subsequent 5 s epoch was selected. The 
epochs were analyzed using a user defined frequency allocation 
into bands of: delta (1-4 Hz), theta (4-8 Hz), alphal (8-10 Hz), 
alpha2 (10-12 Hz), beta (12-25 Hz), and gamma (25-42 Hz). 

Precautions were taken to ensure the gamma band analysis was 
not confounded by electromyographic (EMG) muscle activity or 
eye-saccades artifacts. All conditions were conducted with eyes- 
closed. Each epoch was visually inspected for artifacts and the cut- 
off frequency for gamma was set well below the EMG frequency 
range, which peaks at 70-80 Hz (Lutz et al, 2004). 

EEG SOURCE ANALYSIS 

Based on the scalp-recorded electric potential distribution, 
the exact low resolution brain electromagnetic tomography 
(eLORETA) software (publicly available free academic software at 
http://www.uzh.ch/keyinst/loreta.htm) was used to compute the 
cortical three-dimensional distribution of current source density 
(CSD). The eLORETA method is a discrete, three-dimensional 
(3D) distributed, linear, weighted minimum norm inverse solu- 
tion. The particular weights used in eLORETA endow the tomog- 
raphy with the property of exact localization to test point sources, 
yielding images of current density with exact localization but low 
spatial resolution (neighboring neuronal sources will be highly 
correlated). The description of the method together with the 
proof of its exact zero-error localization property, are described 
in two papers by Pascual-Marqui (2007, 2009). It is important to 
note that eLORETA has no localization bias even in the presence 
of structured noise which constitutes an improvement over the 
previous tomographies of LORETA (Pascual-Marqui et al, 1994) 
and the standardized version sLORETA (Pascual-Marqui, 2002). 
It is important in the context of assessing length of practice related 



differences in anterior and posterior hubs of the DMN that activ- 
ity in these deep structures can be correctly localized with these 
methods (Pizzagalli et al, 2001; Zumsteg et al., 2006). 

Current eLORETA computations were made using a realis- 
tic head model (Fuchs et al., 2002), using the MNI152 template 
(Mazziotta et al, 2001), with the three-dimensional solution 
space restricted to cortical gray matter, as determined by the prob- 
abilistic Talairach atlas (Lancaster et al., 2000). Standard electrode 
positions on the MNI152 scalp were taken from Jurcak et al. 
(2007) and from Oostenveld and Praamstra (2001). The intrac- 
erebral volume is partitioned in 6239 voxels of 5 x 5 x 5 mm 3 
spatial resolution. These eLORETA images represent the electric 
activity at each voxel in Montreal Neurological Institute (MNI) 
space as the exact magnitude of the estimated current density. 
Anatomical labels and/or Brodmann areas are reported using 
MNI space, with correction to Talairach space (Brett et al, 2002). 

The KEY Institute eLORETA software package was used to 
perform these statistical analyses. The methodology used is non- 
parametric. It is based on estimating, via randomized permuta- 
tion testing, the empirical probability distribution for the value 
of the maximum statistic across all voxels under the null hypoth- 
esis. This methodology corrects for multiple testing (i.e. for the 
collection of tests performed for all voxels, and for all discrete 
frequencies). Due to the non-parametric nature of the method, 
its validity need not rely on any assumption of Gaussianity. The 
reader is referred to Nichols and Holmes (2002) for a detailed 
overview of this methodology. 

RESULTS 

SUBJECTIVE RATINGS 

The subjective reports of "meditation depth" are shown in 
Figure 1 . A low level of "meditation depth" was reported in the 
non-meditation (Calculation) condition, with increasing depth 
reported with progression through the meditation conditions. 
The highest rating for meditation depth was for Meditation 
2 (Mantra), although this may include a duration effect. The 
reported levels of ratings were similar between the groups, but 
with slightly higher levels for the student group (SYS) than the 
teacher group (SYT). 



Meditation 'depth' rating 




Meditation 1 (Body-steadiness) 



Meditation z (Mantra) 



FIGURE 1 | Subjective ratings of meditation "depth." 
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COMPARISON OF GROUPS ACROSS CONDITIONS 

We conducted a single test for CSD differences between SYS and 
SYT groups including all frequency bands (delta, theta, alphal, 
alpha2, beta, and gamma). Voxel by voxel the between group log 
F ratio (equivalent to a f-test statistic) was calculated for each 
frequency band considered simultaneously, resulting in 6239 x 6 
comparisons. Following the procedure of Nichols and Holmes 
(2002) implemented in the KEY Institute eLORETA software 
package, we applied 5000 random permutations to calculate the 
distribution of the F-max voxel statistic (the maximum .F-value 
for the full voxel set). The resulting distribution establishes the 
threshold for significance of the F statistic obtained at individ- 
ual voxels in a way which controls for the family-wise error rate 
due to multiple testing. The p-values presented below represent 
the probability of the obtained F statistic at each voxel under 
the null hypothesis while simultaneously correcting for multiple 
testing. 

Tables 1-3 show for each frequency band the region, 
Brodmann area (BA), Tailairach coordinates and the (absolute) 



maximum voxel F statistic or statistical difference between SYT 
and SYS groups and total number of significant voxels (two 
tailed threshold) for Meditation 1 (Body-steadiness), Meditation 
2 (Mantra), and the non-meditation condition (Calculation), 
respectively. 

DIFFERENCES BETWEEN GROUPS ON SOURCE ACTIVATION 
FREQUENCY AND LOCATION IN MEDITATION CONDITIONS 

The two meditation conditions showed a similar pattern, the SYS 
(intermediate group) having voxels with significantly higher CSD 
in the lower frequencies (theta, alphal) and the SYT (advanced 
group) having voxels with significantly higher CSD in the higher 
frequencies (beta, gamma). The SYS group had significantly 
higher CSD voxels in the alpha2 band in Meditation 1 (Body- 
steadiness) but not Meditation 2 (Mantra). No significant voxel 
differences were found in either meditation condition in the delta 
band. 

For the SYS group compared to the SYT group the great- 
est number of voxels with significantly higher CSD values were 



Table 1 | SYT 


> SYS comparison for Meditation 1 


— (Body-steadiness). 








Band 


Max/min region 


Lobe 


Max/min BA 


Tailairach coordinates 


Max/min F-value 


No. of voxels p < 0.05 


Delta 
Theta 


Inf parietal lobule 
Inf parietal lobule 


Parietal 
Parietal 


R40 
R40 


40, -40, 60 
45, -45, 60 


-0.729 ns 
-1.32** 


490 


Alphal 


Precentral gyrus 


Frontal 




55, -10, 45 


-2.01** 




Alpha2 

Beta 


Inf parietal lobule 
Sub— gyral 


Parietal 
Temporal 


R1 
R20 


45, -30, 65 
40, -10, -25 


-1.08* 
1.03* 


30 
33 


Gamma 


Fusiform gyrus 


Temporal 


R20 


40, -15, -30 


1.65** 







»p < 0.05, *»p < 0.01. 



Table 2 | SYT > SYS comparison for Meditation 2— (Mantra). 



Band 


Max/min region 


Lobe 


Max/min BA 


Tailairach coordinates 


Max/min F-value 


No. of voxels p < 0.05 


Delta 
Theta 


Ant cingulate 
Inf parietal lobule 


Limbic 
Parietal 


L32 
R40 


-10, 35, -5 
50, -50, 55 


0.822 ns 
-1.52** 


613 


Alphal 


Precentral gyrus 


Frontal 


R4 


50, -10, 50 


-1.97** 




Alpha2 
Beta 


Sup temp gyrus 
Rectal gyrus 


Temporal 
Frontal 


R22 
L11 


45, -20, 0 
-10, 40, -25 


1.10 ns 
1.30* 


917 


Gamma 


Insula 


Sub-lobar 


R13 


30, 20, 15 


1.92** 







*p < 0.05, **p < 0.01. 



Table 3 | SYT > SYS comparison for control condition (Mental calculation). 



Band 


Max/min region 


Lobe 


Max/min BA 


Tailairach coordinates 


Max/min F-value 


No. of voxels p < 0.05 


Delta 


Sup temp gyrus 


Temporal 


R41 


55, -25, 5 


-1.72** 


390 


Theta 


Sup temp gyrus 


Temporal 


R22 


65, -20, 0 


-1.07 ns 




Alphal 


Precentral gyrus 


Frontal 


R6 


45, -10, 40 


-1.96** 




Alpha2 


Inf parietal lobule 


Parietal 


R40 


50, -50, 55 


-1.43* 


123 


Beta 


Sub— gyral 


Temporal 


L20 


-45, -10, -25 


1.07 ns 




Gamma 


Fusiform gyrus 


Temporal 


R20 


40, -15, -30 


1.52* 







*p < 0.05, **p < 0.01. 

Positive value for F indicates SYT shows higher CSD than STS, negative value for F indicates SYS has higher CSD than SYT 
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in the alpha 1 band, predominantly in the right hemisphere. For 
both meditation practices, the three Brodmann areas in which 
the highest alpha 1 voxel F- values were located were, in rank 
order of that statistic, right BA4 (precentral gyrus), BA6 (pre- 
central gyrus), and BA3 (primary somatosensory cortex). (See 
Figures 2, 3). 

For the SYT group compared to the SYS group, the great- 
est number of voxels with significantly higher CSD occurred in 
the gamma band, predominantly in the right hemisphere. In 
Meditation 1 (Body-steadiness), the highest voxel statistic val- 
ues were located in right BA20 (fusiform gyrus), BA21 (middle 
temporal gyrus), BA13 (insula), BA38 (the anterior pole of the 
temporal lobe). In Meditation 2 (Mantra), the highest voxel statis- 
tic values were located in right BA13 (insula), BA45, and BA47 
(inferior prefrontal gyrus). (See Figures 5, 6). 

DIFFERENCES BETWEEN GROUPS ON SOURCE ACTIVATION 
FREQUENCY AND LOCATION IN CALCULATION CONDITION 

Significant differences were found for voxels across the delta, 
alpha 1, alpha2, and gamma frequency bands, with the SYS group 
showing significantly higher CSD voxels than the SYT group in 
delta, alphal, and alpha2 bands. In these low frequencies, the 
locations of significant voxels were entirely in the right hemi- 
sphere, the voxels with the highest F-values being located in 
rank order in right BA41, 22, 21 in delta; in right BA6, 4, 3 in 



alphal; and in right BA40, 7, 19 in alpha2. (See Figure 4 for 
alphal). 

In the high frequency gamma band, the SYT group showed 
voxels with significantly more CSD than the SYS group. Once 
again these voxels were found entirely in the right hemisphere, 
the voxels with the highest F-values being located in rank order in 
right BA20, 13, 21. (See Figure 7 for gamma). 

CORTICAL SOURCES OF DIFFERENCES BETWEEN GROUPS IN HIGH 
AND LOW FREQUENCY BANDS 

Across all conditions the number of significant differences was 
consistently greatest in the gamma band for the high fre- 
quencies and the alphal band for the low frequencies (high- 
lighted in Tables 1-3 above). In order to examine the spe- 
cific regional pattern of these high and low frequency group 
differences, the number of significant voxels in each (right 
and left) Brodmann Area are presented for each experi- 
mental condition in the alphal band and gamma bands in 
Tables 4, 5. 

DISCUSSION 
HYPOTHESES 1 AND 2 

The first prediction that during Yoga meditation student (trainee- 
teacher) meditators would show greater CSD in the low fre- 
quency bands was supported in the theta, alphal and alpha2 




FIGURE 2 | Differences in CSD between groups in alphal band— Meditation 1 (Body-steadiness) Increased alphal activity (blue) in SYS compared to 
SYT in right BA4 (precentral gyrus). 
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FIGURE 3 | Differences in CSD between groups in alphal band— Meditation 2 (Mantra) Increased alphal activity (blue) in SYS compared to SYT in 
right BA4 (precentral gyrus). 
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FIGURE 4 | Differences in CSD between groups in alphal band— Non-meditation (Calculation) Increased alphal activity (blue) in SYS compared to 
SYT in right BA6 (precentral gyrus). 
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FIGURE 5 | Differences in CSD between groups in gamma band in Meditation 1 (Body-steadiness) Increased gamma activity (yellow) in SYT 
compared to SYS in right BA20 (fusiform gyrus). 
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FIGURE 6 | Differences in CSD between groups in gamma band in Meditation 2 (Mantra) Increased gamma activity (yellow) in SYT compared to SYS 
in right BA13 (insula). 



bands for Meditation 1 (Body-steadiness) and theta and alpha 
1 for Meditation 2 (Mantra). In both meditations the number 
of significant voxels at the lower frequencies was overwhelmingly 
greatest in the alphal band (See Tables 1, 2). 

The second prediction that advanced practitioners would show 
greater CSD in the high frequency bands was supported for 
beta and gamma frequency bands in both meditation conditions. 
For the higher frequencies the number of significant voxels was 



overwhelmingly greatest in the gamma band (See Tables 1, 2). 
This is the first study to show enhanced gamma band activation 
in advanced Western meditators practicing in the Yoga tradition, 
compared to less experienced practitioners. 

More broadly the P-max voxel values for SYT vs. SYS were 
negative for all low frequencies (delta to alpha2) in all experi- 
mental conditions (including mental calculation) and positive for 
all higher frequency bands (beta and gamma), (See Tables 1-3). 
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FIGURE 7 | Differences in CSD between groups in gamma band in Non-meditation (Calculation) Increased gamma activity (yellow) in SYT compared 
to SYS in right BA20 (fusiform gyrus). 



It appears then that patterns of increased (low frequency) EEG 
activity in Yoga meditation students are consistent with the 
tradition based expectations of pratyahara (somatosensory with- 
drawal), evidenced by significantly greater low frequency activ- 
ity corresponding to inhibition of cortical processing, in those 
regions mediating external sensory and motor processing. By 
contrast, patterns of increased (high frequency) EEG activity in 
advanced Yoga teachers most likely correspond to activity in brain 
regions recruited in the conscious states of concentration and 
absorption, which emerge in advanced practitioners. 

While the relationship of the neural networks engaged by 
meditation practice and hypnotic induction remains to be deter- 
mined it is intriguing to note that a recent study of the neu- 
rophenomenology of neutral hypnosis (what some would term 
trance) found that for those with high hypnotic susceptibility self- 
reported hypnotic depth correlated significantly with EEG activ- 
ity measures uniquely in beta and gamma band frequency ranges 
(Cardena et al., 2013). At the same time the high (but not low 
or medium) susceptible group reported spontaneous exceptional 
experiences of positive affect and/or self-transcendence. 

HYPOTHESES 3 AND 4 

Tables 4, 5 show that for the alphal and gamma bands, respec- 
tively, only right hemisphere Brodmann Areas show significant 
differences between SYT and SYS in the non-meditation (count- 
ing backwards by four) condition. Looking across the rows 
of these tables these same regions are also found to have the 
greatest number of significant voxels in the "Body-steadiness" 
and "Mantra" meditation conditions (typically with a greatly 
increased number of significant voxels). We interpret this pattern 
as support for our third prediction of trait differences in baseline 
neural activity brought about by long term meditation practice 
in those regions most engaged by this practice. We further inter- 
pret these findings as evidence for a core right-sided (rather than 
midline) network (or networks if the alphal and gamma results 
are considered separately) that is progressively modulated over 
the course of Yoga meditation practice. The predominantly right 
lateralized location of the cortical sources where activation dif- 
ferentiates between the STS and SYT groups can be observed for 
both alphal and gamma in each of the experimental conditions 
in Figures 2, 3, 5, 6. Thus, in the case of Yoga meditation practice 



(as distinct from recent studies of mindfulness/Buddhist medita- 
tion practices) the midline nodes of the DMN, although included 
in regions of significant voxel differences, do not appear to be 
a principal locus of practice-related changes in cortical activity 
(contradicting the fourth prediction above). Our result also dif- 
fers from the CSD changes in midline DMN structures recently 
reported in "concentrative" (breath-focused) meditation with 
"intermediate level" (mean 4 years) practitioners of unspecified 
tradition compared to controls (Lavallee et al, 2011). It also sug- 
gests an important difference (in advanced meditators) with the 
neutral hypnotic state in so much as hypnotic induction has been 
found to reduce (fMRI measured) anterior DMN activity in high 
hypnotically susceptible participants (McGeown et al., 2009). 

REGIONAL DIFFERENCES IN EXPERIMENTAL CONDITIONS RELATED TO 
YOGA MEDITATION PROFICIENCY ALPHA1 BAND 

In the alphal band, significant voxel differences always showed 
greater activity in the student than the teacher meditators. In 
all conditions, the most significant voxel F-values were always 
located in the right somatosensory (BA1, 2, 3, 5), motor (BA4) 
and premotor cortex (BA6). In the mental calculation condition, 
a core network of right sided regions can be identified. This was 
comprised, in addition to those regions identified immediately 
above, of right occipital (BA17, 18, 19), right parietal (BA7, 39, 
40) and right superior frontal gyrus (BA8, 9). These same regions 
remain the core of the most significant proficiency related dif- 
ferences in both meditation conditions. The key changes are an 
increase in the number of significant voxels recruited in each 
region and a spread of significant voxels to the homologous left 
sided Brodmann Areas, while always retaining a strong imbal- 
ance toward right sided voxels (see Table 4). Following current 
models of the functional role of task related topographic dif- 
ferences in alpha band activity, these effects may be interpreted 
as trained task specific patterns of functional cortical inhibi- 
tion (Klimesch et al., 2007). Such an interpretation of present 
alphal findings fits remarkably closely with the specific effects of 
the inhibitory thalamo-cortical mechanism proposed by Austin 
(2013) to account for states of deep absorption corresponding 
to Yoga meditation practices. However, in Austin's model these 
effects would be expected to be greatest in advanced rather than 
intermediate practitioners. 
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Table 4 | Number of significant voxels in each Brodmann Area for 
each condition in alphal band. 



No. signif voxels 
BA 


Alphal band (8-10 Hz) 


SYS 


> SYT 


Calculation 


Body 


Mantra 


L R L 


R 


L 


R 


I 


9 


9 




Q 


Z 


47 


49 


Z. 


ro 
Do 


Q 
O 


59 12 


72 


17 
1 / 


/ u 


A 

4 


54 12 


63 


1 Q 
I o 


R 1 
D I 


O 


12 27 


50 






c 
D 


170 38 


265 


04 


0K7 
ZD / 


/ 


130 92 


208 


1 1 Q 
I I O 


OnQ 
ZUo 


Q 
O 


33 


64 


1 


7R 
/O 


Q 

y 


55 


84 




oy 


17 
I / 


8 


40 




1 O 
I o 


I O 


25 


158 




40 


i y 


64 6 


145 


2 


98 


Z J 


3 


6 




Z 


Z4 


1 


34 


14 


44 


30 


1 


18 




6 


o I 


12 12 


90 


27 


69 


32 




12 


2 


20 


37 




14 




3 


39 


28 


58 




55 


40 


143 


163 


3 


149 


43 




3 






44 




3 






45 


1 


6 




1 


46 




4 






Total 


0 856 199 


1627 


310 


1371 



Brodmann areas with no voxels omitted from table. 
GAMMA BAND 

In the gamma band, significant voxel differences always showed 
greater activity in the teacher than the student meditators. In 
the non-meditation (mental calculation) condition all significant 
voxel differences were found in the right anterior temporal lobe 
(principally BA20, 21, 22, 38) and right ventral prefrontal cortex 
(principally BA13 and BA44). All these regions showed a great 
expansion of significant voxel counts in the body-steadiness and 
mantra meditations. 

In the "Body-steadiness" meditation, 1704 out of a possible 
6239 cortical voxels showed significant proficiency differences and 
in the "Mantra" meditation this rose to 2090 significant voxel dif- 
ferences. A right sided bias continued to be observed for both 
meditations but was much less extreme than in the alphal band. 
Additional regions showed significant greater CSD in gamma 
for the long term meditators in both meditations (principally in 
BA10, 11,45,47) (see Table 5). 

This result adds to the converging lines of enquiry regarding 
the role of meditation-specific increases in gamma band activ- 
ity in "advanced" practitioners, adding evidence from Western 
meditators practicing in an integral Yoga spiritual tradition. The 



Table 5 | Number of significant voxels in each Brodmann Area for 
each condition in gamma band. 
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Brodmann areas with no voxels omitted from table. 

teacher group in this study (mean experience 30 years) displayed 
striking increases in gamma band activity similar to studies of 
"advanced" Tibetan Buddhist meditators (Lutz et al., 2004) and 
Western Buddhist Vipassana meditators with mean experience 20 
years (Cahn et al., 2010). Although those studies did not esti- 
mate the activity of cortical sources, the high level of gamma 
activity during meditation in advanced practitioners requires 
widespread synchronization throughout extensive cortico-cortico 
and cortico-thalamic neural networks (Lutz et al., 2004). 

In Meditation 1, the "Body-steadiness" meditation, the most 
significant voxel differences were located in the right ante- 
rior temporal lobe and the insula. This result aligns with that 
of an advanced Tibetan Buddhist meditator showing increased 
gamma activity in right mid temporal gyrus (BA21) in the 
comparisons of sensory-focused "visualization" meditation vs. 
the verbal-focussed "mantra" meditation and also in a "self- 
reconstruction" vs. "self-dissolution" meditation (Lehmann et al., 
2001). In Meditation 2, the "Mantra" meditation, the most signif- 
icant sources were located in the right insula and right inferior 
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frontal gyrus. The right insula has been linked with a more 
detached and objective awareness of interoceptive sensory events 
(Farb et al., 2007) involved in the shift from "narrative" to "expe- 
riential" self- awareness. Deen et al. (2011) identify three distinct 
clusters of functional connectivity with network hubs in pos- 
terior, mid, and anterior insula, respectively, suggesting that a 
further parcelation of anatomical and function subregions within 
the insula will be required to fully understand the role (or roles) 
it plays in Yoga and Buddhist tradition meditative states. It may 
be of note that a major white fiber tract, the uncinated fasciculis, 
mediates bidirectional connections between the high level associ- 
ation cortex of the anterior temporal lobe and the inferior frontal 
gyrus (Von der Heide et al., 2013), the only such structure which 
continues to develop beyond the age of 30 (Lebel et al., 2008). 

STRENGTHS OF THIS STUDY 

A strength of this study is the advanced level of meditation expe- 
rience in a sample of Western Yoga practitioners following an 
identical spiritual tradition. All of the "advanced" group were san- 
nyasin (initiated) disciples of Swami Satyananda who had spent 
considerable time in ashrams, both in Australia and India. They 
were all Yoga teachers. Studies reporting results from large groups 
of advanced meditation practitioners often combine meditators 
from different traditions into a single experimental group. Such 
differences are not considered equivalent by the members of these 
traditions themselves and any such equivalence must be demon- 
strated rather than assumed. We have chosen not to make this 
assumption even at the cost of a lower sample size in order to 
maintain the validity of our analysis. A further strength was the 
attention given to preserving the ecological validity of the medi- 
tation states attained, the study being conducted in a meditation 
setting conducive to the participants' usual practice. The intru- 
sion of equipment into this setting was minimized, with the 
EEG cap and small radio transmitter being the only items in the 
meditation room. 

As with all EEG source analysis methods, eLORETA results 
are highly susceptible to distortion by artifacts (Pizzagalli, 2007), 
therefore the quality of data preprocessing is essential to the 
validity of the results obtained. We adopted a recently devel- 
oped (Nolan et al., 2010) ICA component rejection procedure for 
artifact cleaning developed by the Neural Engineering Group at 
Trinity College Dublin. This is one of a new generation of ICA 
artifact correction procedures which is guided by a series of objec- 
tive component parameters rather than the momentary choice 
of the individual experimenter. We believe the effective utiliza- 
tion of this method for preparing large volumes of clean data was 
an essential factor contributing to the validity of this research. 
We strongly recommend the future adoption of this or similar 
method by related research programs. 

LIMITATIONS OF THIS STUDY 

Limitations of the study include the small sample size (although 
that did not prevent significant findings) and the potential 
confound due to significant age difference between the sub- 
ject groups. A larger sample size is desirable to enhance the 
power and generalizability of the current results and would allow 
the additional possibility of parametric testing (Thatcher et al., 



2005) . However, the non-parametric statistical mapping method 
adopted here is well suited to the limitations of the current design 
and sample size. In fact, one of the worked examples offered by 
Nichols and Holmes (2002, p. 15) is a PET study with 12 sub- 
jects equally divided amongst 2 testing conditions which directly 
reflects the constraints of the present study. In addition the ran- 
domization distribution is known to be overly conservative when 
the effect is distributed, as occurs in the present results (Troendle, 
1995). Therefore, we believe that there is a strong justification for 
reporting the group effects identified here. 

EEG spectral power density changes dramatically with age 
until early adulthood. From that point changes are slow in healthy 
aging. Power in lower frequency bands (delta through alpha) 
decreases throughout this time period (Rossini et al., 2007). The 
age range of our subjects extends from 30 to 63, so it could be that 
significantly greater activity in the younger group in these lower 
bands is due to normal aging rather than years of meditation 
practice. 

A comparison of the effects of normal aging to the presently 
observed differences between the SYS (younger) and SYT (older) 
studies of EEG spectral source changes across this part of the lifes- 
pan shows that delta decreases in occipital sources (Babiloni et al., 

2006) , but we only found significant delta differences in one con- 
dition (counting) and this was maximal in right BA41 (primary 
auditory cortex). In normal aging alpha 1 and alpha2 decreases in 
occipital, parietal, temporal and limbic cortices (Babiloni et al., 
2006). We found significantly lower alphal in occipital and pari- 
etal cortex but not in temporal or limbic cortex. Alpha also 
increases with age in frontal cortex (Basar, 2012) but we found 
significant alpha decreases in superior and middle frontal gyri in 
the older group and no significant differences in inferior frontal 
gyri. Once again, the pattern of our findings does not fit with what 
would be expected if they were due to normal aging. 

Despite available data there is no clear evidence for changes 
due to normal aging through middle age in the high frequency 
bands. Therefore, we do not consider age differences to be a plau- 
sible alternative explanation for those high frequency effects. A 
confound due to the age difference between the groups fails to 
account for why the set of BA containing significantly greater 
low frequency voxel activation in the SYS group forms the com- 
plement to the set of BA containing significantly greater high 
frequency voxel activation in the SYT group. Nor does it explain 
why a core pattern found in the non-meditation Calculation 
condition is maintained but greatly expanded in the Body Scan 
and Mantra meditation conditions. Lutz et al. (2004), in a study 
with advanced and novice meditators with age means of 49 years 
and 2 1 years, found hours of practice, but not age, significantly 
predicted gamma in their baseline condition. 

The additional analysis of age matched subgroups goes some 
way toward addressing this issue with significant F-max voxel 
findings in both alphal and gamma bands located in the same 
regions as the original group analysis. In this case the further 
reduced sample size points to the robustness of these effects. 
Notwithstanding these considerations practice related differences 
in age will be intrinsic to group membership whenever truly 
"advanced" (such as the 30 year group in the present study) 
are compared with medium term practice groups. Now that 
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specific regions of interest have been identified parametric meth- 
ods which enable statistical control of age effects present a viable 
option to address this issue in future (larger n) studies. 

The density (number) of electrodes in the recording array 
also places limitations on the resolution of source localization 
(Pizzagalli, 2007). Dense array recordings are so far limited to lab- 
oratory settings, unlike the present study. However, successful val- 
idation studies of LORETA, sLORETA and eLORETA use only the 
original 19 electrodes of the 10-20 system (Pascual-Marqui et al., 
2002). Particularly for non-evoked recording conditions, 19 chan- 
nel recordings are considered adequate and are widely reported 
for the LORETA family of source analysis methods (Babiloni 
et al, 2014). By comparison, the current study employed 25 active 
recording electrodes. A further possible limitation was the fixed 
sequence of conditions rather than a systematically counterbal- 
anced design. For investigation of deeper meditation states, the 
benefits of experimental manipulation of conditions and intrusive 
measuring instruments must be balanced by the need to obtain an 
authentic meditation experience. 

CONCLUSION 

This study is the first to report enhanced gamma band activation 
in advanced Western meditators practicing in the (Satyananda) 
Yoga tradition, compared to less experienced (intermediate) prac- 
titioners. It found significant differences in EEG frequency band 
sources between "experienced" (mean 4 years) and "advanced" 
(mean 30 years) meditators in two meditation conditions (Body- 
steadiness and Mantra) and a mental calculation condition. 
The findings, strongest in alphal and gamma bands, consis- 
tently support the frequency band hypothesis proposed by Fell 
et al. (2010) of enhanced low and high frequency band effects, 
respectively, for intermediate and advanced levels of medita- 
tion experience. It further adds a finding of increased gamma 
band activity in advanced western Yoga meditators to similar 
findings in advanced Buddhist meditators. In addition, the Body- 
steadiness meditation condition in the present study appears to 
have much in common with Buddhist mindfulness of breath- 
ing practice. Both practices direct awareness to the felt quality of 
bodily sensations. The present findings then call for much closer 
scrutiny of the distinction, at least in the long term, between 
the functional networks trained by concentrative and mind- 
fulness meditation methods. Those regions which consistently 
differentiated between the groups in the low and high frequency 
bands across all meditation conditions showed similar (but 
reduced) differences in the non-meditation (counting backwards) 
condition. 

This study extends previous EEG findings by estimating cor- 
tical gray matter sources for frequency band specific meditation 
training effects. The cortical loci of high and low frequency 
band meditation training effects are both specific and distinct. 
Approximately one third of available cortical gray matter vox- 
els showed significantly greater gamma band activity in mantra 
meditation in advanced (SYT) than intermediate (SYS) med- 
itators. Despite this, there is little overlap with those regions 
which show increased alphal in the student meditators. In fact 
the Brodmann areas in which significant voxel differences are 
found during meditation in the alphal (suggesting functional 



inhibition) and gamma bands (suggesting integration into con- 
sciousness) form almost absolute complement subsets of the full 
set of cytoarchitecturally defined gray matter regions. We propose 
that selective inhibition of a right lateralized network comprising 
visual, somatosensory and body-world self-representations cor- 
responds to the earlier stages of sensory withdrawal and stripping 
away of the "outer" onion-like layers of self as described in tra- 
ditional Yoga literature. The subsequent emergence of conscious 
states specific to advanced practitioners requires both the dis- 
engagement from these self-world representational systems and 
the development of widespread gamma synchronization through- 
out anterior temporal and ventral prefrontal cortical regions 
extending from a right sided core network incorporating anterior 
temporal lobe and insula. 
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